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We explore static spherically symmetric stars in the Gauss-Bonnet gravity with- 
out cosmological constant, and present an exact internal solution which attaches to 
the exterior vacuum solution outside stars. It turns out that the presence of the 
Gauss-Bonnet term with a positive coupling constant completely changes thermal 
and gravitational energies, and the upper bound of red shift of spectral lines from 
the surface of stars. Unlike in general relativity, the upper bound of red shift is 
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^j^- dependent on the density of stars in our case. Moreover, we have proven that two 

^f^ ', theorems for judging the stability of equilibrium of stars in general relativity can be 

> 

^\\ ' hold in Gauss-Bonnet gravity. 
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I. INTRODUCTION 

Up to now, many quantum theories of gravity have been proposed. The most promising 
candidate among them is the superstring/M-theory. In string theory, extra dimensions were 



H ' promoted from an interesting curiosity to a theoretical necessity since superstring theory 

requires an eleven-dimensional space-time to be consistent from the quantum point of view 



B 



^. It caused a renewed interest about the general relativity in more than 4 dimensions. 
On the other hand, in recent years another renewed interest has grown in higher order 
gravity, which involves higher derivative curvature terms. Among the higher curvature 
gravities, the most extensively studied theory is so-called Gauss-Bonnet gravity |5|, |6| , which 
can be naturally emerged when we want to generalize Einsteins theory in higher dimension 
by keeping all characteristics of usual general relativity excepting the linear dependence of 
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Riemann tensor. In addition, Gauss-Bonnet term occurs in the effective low-energy action of 



superstring theory 7H10l|. Hence, the Gauss-Bonnet gravity provides a promising framework 
to study curvature corrections to the Einstein-Hilbert action. In Gauss-Bonnet gravity, the 
exact static and spherically symmetric black hole solutions were found in Ill4l3l|. their 
thermodynamics have been investigated 14J-ll7l| and the slowly rotating black hole solutions 



have been obtained in 



18 



2l|. 



As well known, in general relativity, there is an exact solution of gravitational fields 
with a spherically symmetric static star, which takes the form of the exterior Schwarzschild 



solution of vacuum outside the star 



22| . It implies that the exterior Schwarzschild solution 



can be naturally emerged by a suitable distribution of matter. Thus, we are convinced of the 
existence of the Schwarzschild black hole, since it is formed from the gravitational collapse 
of a heavy star. An interesting question will be, what's the solution of gravitational fields 
with a spherically symmetric static star in higher dimensional Gauss-Bonnet gravity, which 
relates to the solution of black holes, if matter can disperse to the extra dimension? 

The objective of this paper is to find a solution and discuss the physical properties of 
static spherically symmetric stars in Gauss-Bonnet gravity without cosmological constant. 
The dimension of space-time is assumed as D > 5 since Gauss-Bonnet term yields nontrivial 
dynamics. We try to investigate the spherically symmetric system consisting of static perfect 
fluid, and give an exact solution. In order to find the effect of Gauss-Bonnet term, we discuss 
thermal and gravitational energies, the distribution of pressure of stars, and the upper bound 
of red shift of spectral lines from the surface of them. Meanwhile, we also study that whether 
two theorems for judging the stability of equilibrium in general relativity can be hold in 
D > 5 Gauss- Bonnet gravity. 

This paper is organized as follows. In Section 2, we give the equations for stars' structure 
in D > 5 Gauss-Bonnet gravity without the cosmological constant, and find a solution. In 
Section 3, we investigate the effects of Gauss-Bonnet term. The stability of equilibrium of 
stars is demonstrated in Section 4. Some remarks are given in Section 5. 



II. DIFFERENTIAL EQUATIONS AND SOLUTIONS FOR STAR STRUCTURE 

The action of Gauss-Bonnet gravity in D- dimensional space-time can be written in the 
form of 

•S* = -— / d^Xy/^{R + aC2) + Smatter, (l) 

where the couphng constant a can be regarded as the inverse string tension and we assume 
a > in this paper. Here the second order (Gauss- Bonnet) term of Lagrangian is given by 

Varying the action Eq. ([T]), we obtain the field equations as 

Gab = GI;^ + «G1? = KTab, (2) 

where 

Cab ^ ^ab — -Rgab, 

^ab ~ '^{~RacdeR ^'^b ~ '^RacbdR'^ — '^RacR'^b + RRab) — ■::^2gab- 

For future simplicity, we take coefficient a = (£)_3)%_4) • 

The metric for a static spherically symmetric system will be taken in the form of 

cis^ = -B{r)dt^ + A{r)dr'^ + r^dnl_^, (3) 

where A and B are functions of r, and dQ'jj_2 represents the line element oi a {D — 2)- 
dimensional unit sphere with f2£)_2 = 2n'^^^^^^'^/T[(D — l)/2]. The energy-momentum tensor 
is assumed to be that for a perfect fluid 

T^,u = (p + p)U^U^ + pg^,u, (4) 

with p the proper pressure, p the proper total energy density, and U^ the velocity four- vector, 
defined so that UJJ^ = —1. Since the fiuid is at rest, we take If^ = J^ — Si". We assume p 
and p are positive throughout the star. Our assumptions of time independence and spherical 
symmetry imply that p and p are functions only of the radial coordinate r. 



Plugging Eq. (|3]) and Eq. @ into Eq. ([2]), we find that the equations read 

, _ D-2 a[{D -5){A- if A + 2{A - l)A'r] + {D - 3){A - 1)AV + AA'r^ 

= K,p, (5) 

_ D - 2 a[-{D -5){A- ifB + 2{A - l)B'r] - {D - 3){A - l)ABr^ + AB'r^ 

= f^P, (6) 

where an prime denote the derivative with respect to r. Here we omit the third equation. 
Instead, we use the equation for hydrostatic equihbrium which is equivalent to the equation 
of energy-momentum conservation 

B' _ 2p' 
B P + p' 

From Eq. Q we find that the solution with A{0) finite is 



(7) 



Mr) = , / „.. =. (8) 
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where the mass function M{r) is defined as M{r) = i^D-2 C p{r')r'^^^'^^dr'. Here we have 
abandoned another solution since it doesn't recover the solution in general relativity in the 
limit a — > 0. The function M{r) isn't the total energy of the matter since it is defined as 
the integral of the energy density p{r) of matter alone thus doesn't include the energy of 
the gravitational field. The thermal and gravitational energies of the star will be discussed 
in the next section. Using Eq. IQ, Eq. ([7]) and Eq. ([8]), we obtain 






2kM I 2(D-5)kM 



p=-r(p + p) , , ■, (9) 

A?. 1 I 2kM m .f__ / r^ I 2nM N 
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and the solution with Bioo) = 1 is 



,_ ^r)-U9 /J- J 2kM I 2(D-5)kM 1 



D-2^ ' ^y 4a^ ' (D-2)5ni5,2rO-i ' (D-2)nB_2rO-i " , , ,^^, 

ij(r ) = exp — / r nr\ (111) 



y 4a2 ^ (D-2)aQ.r,~2r'^-^^ 2a \l icP- ' {D-2)aQ.D-2r^'-' ' 

Outside the star, p(r) and p{r) vanish, and M{r) is the constant M{K), therefore Eq. ([8 
and Eq. flTOj) give 



BM^^-V)^i + ^-J^+(jj3^|^^ (11) 



for arbitrarily r > R. The Eq. (ITT]) is the exterior solution of the static spherically symmetric 
black hole in fiat space. Hence, if matter can disperse to the extra dimension, in the D > 5 
spherically symmetric space-time, there is an interior solution which describes a perfect 
fluid relating to the exterior solution of the static spherically symmetric black hole. Thus, 
Eq. ( ITTl) makes the solution of black holes more believable since such exterior solution can 
be naturally emerged by a suitable distribution of matter. 

The Eq. (|9]) is the generalization of Tolman-Oppenheimer-Volkoff (TOV) equation in the 
D > 5 EGB theory. One can check that in the general relativity limit a — )■ 0, such an 
equation in 5D space-time will be written in the form of 

_l_ {D-3)kM 

p' = -''(p + p) j'1/°-Cm ■ (12) 



which is similar to the TOV equation in 4D case |24J-l26| . As in general relativity, we can use 
this equation to obtain p{r), M{r), p{r) throughout a static spherically symmetric star in 
the 5D space-time. The pressure p may in general be expressed as a function of the density 
p, the entropy per nucleon s, and the chemical composition. We can assume that the entropy 
per nucleon s does not vary throughout the star, and the star we consider have a chemical 
composition that is constant throughout for simplicity, that is, p{r) may be regarded as 
a function of p{r) alone. The definition of M{r) provides an initial condition M(0) = 0. 
Eq. ([9]) together with an equation of state giving p{p), serve to determine p{r), M{r), p{r) 
of the star for arbitrarily r, once we specify the other initial condition that the value of p(0). 
Eq. (J9]) should be integrated out from the center of the star, until p{p{r)) drops to zero at 
some point r = R, which we then interpret as the radius of the star. 

III. EFFECT OF GAUSS-BONNET TERM 



This section is devoted to the influence on a spherically symmetric star caused by Gauss- 
Bonnet term. First we investigate thermal and gravitational energies of the star. As in 
general relativity, we can compare M{R) with the energy Mq that the matter of the star 
would have dispersed to infinity in D > 5 space-time. This is simply Mq = mj^N , where 



rrii^ is the rest-mass of a nucleoli, A^ is the number of nucleons in the star which is given by 
N= [ ^f^J%drdOi ■ ■ ■ dOD-2 = ^D-2 I ^A{r)B{r)J%r''-^dr, 

JT. Jo 

where J^ is the conserved nucleon current, 6*1, ■• ■ , 9d-2 are spherical coordinates on Sd-2, 
S is a Cauchy surface of the spacetime. J^ can be expressed in terms of the proper nucleon 

7/' _ , /r 70 



number density n, which is n = —U^J'^ = VBJ^. The expression of A^ then becomes 



N = nD-2 / ^/A(^n{r)r^-^dr. (13) 

Jo 

The proper number density n{r) and A^ are fixed for a star with a given constant s and 
chemical composition, once we choose p(0). The internal energy of the star is now defined by 
E = M — itlnN . We can define the proper internal material energy density as e(r) = p{r) — 
mNn{r), and decompose the internal energy E as the sum of the thermal and gravitational 
energies T and V, respectively 

rR 

T = nD-2 / VAe{r)r^-'^dr, 
Jo 

j-R 

V = nD-2 / (1 - V^)pir)r^'^dr. (14) 

Jo 



Using Eq. (|8]), we obtain that 
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V = nD-2 I (1 - ^ )p{r)r'''Mr. (15) 

/l I Ili _ / r4 I 2kM 
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It is straightforward to check that 

dA \ (, 1 



da. r,~ I r'^ , 2kM 



(l ==) <0, (16) 

V 1 I r2 / r* I 2kM J 



OA, / »■ _| -^«™ \ 1 _|_ £1 _ / '' J ^_ 

^"Y 45^ ^ (D-2)aOD_2r"-5 "^ "•" 2a Y ^^^ (D-2)ar2D_2r"-5 

which leads to dT /da < and dV/da > 0. Hence, Gauss-Bonnet term with positive a will 
decreases the thermal energy and increases the gravitational energy. 

Next, we consider a stars with uniform density, consisting of incompressible fluids with 
p = canst. In this case Eq. (Q may be written as 



2 . r, 1 I 2kp . 2{D-5)kp 1 

D-2P^^y 4a2 ^ (D-l)(D-2)a^ (D-l)(D-2) a 

p (r) = -r(p + p) ^ ^^^ -. — / 1 2k^ ^' ^ ^ 

^'^V 432 + (D-l)(D-2)a\-^ + ^ (23 ~ V 43^ + (D-l)iD-2)a>\ 



The solution with p{R) = is 

here we define 

$ = {D-2). 



^ P£_P4HM, (18) 



2/tp {D - 5)/tp D-2 



4^2 {D-l){D-2)dt D-1 2a 



1 / 1 , 2«:p 



2(5 V^a^ (L>-l)(D-2)a' 



/i + ei?2 

Such solution doesn't make sense for all values of p and R. The pressure will become infinite 
at a point Too where 

2 _ (1 + eR^)K^p^ - $2 

'^oo ~ 0$2 • 

An infinity in p{r) can't be blamed on an injudicious choice of coordinate system since the 
pressure is a scalar. In order to ensure that p{r) isn't singular throughout the star, we find 
that 

i + eR'> ^, (19) 



which leads to r^ < 0. 

Hence, the function B{R) = l + Oi?^ has a lower limit. This result has profound infiuence 
on the red shift of spectral lines from the surface of stars. The red shift z will have an upper 
bound Zmax since z = ^ = B'^i^R) — 1, which gives 

Zma. = ^ - 1- (20) 

In the general relativistic limit, we have \m\ Zmax = 2/(-D — 3). Obviously, non-negative 

a— !>0 

coefficient a will change the upper bound of the red shift, which is given under in FIG.l. 
The most interesting consequence of the of the second order curvature corrections is that 



Zmax is dependent on the value of p while it is a constant in general relativity case 27 1. 

IV. STABILITY OF EQUILIBRIUM 
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{a)D = 5 (h)D = 6 

FIG. 1: Relation between a and the upper bound of z^ax in 5D and 6D Gauss-Bonnet gravity: 
the sohd and dotted hnes representing p = 2 and p = 1 respectively. In general relativity limit 
a —7- two lines have a same limit.. 

Our solution represents an equilibrium state of the star, but it may be stable or unstable. 
For most purposes, we only concerned with the stable solution. In order to investigate the 
stability of a particular configuration, it would be necessary to compute the frequencies a;„ of 
all normal modes of the configuration. If the frequency Un has a positive imaginary part, the 
factor exp{—iUnt) would grow exponentially, and the system would be unstable. However, 
it is often possible to obtain from the equilibrium solution alone whether the corresponding 
configuration is stable, since we have following two theorems: 

1. A star consisting of a perfect fluid with constant chemical composition and entropy 
per nucleon can only pass from stability to instability with respect to some particular radial 
normal mode, at a value of the central density p(0) for which the equilibrium energy E and 
nucleon number N are stationary that 

dE{p{0),s--) 



dp{0) 
dN{p{0),s 
dp{0) 



0, 
0, 



(21) 



by a radial normal mode is meant a mode of oscillation in which the density perturbation 
6p is a function of t and r alone, and in which nuclear reactions, viscosity, heat conduction, 
and radiative energy transfer play no role. 

2. A particular star configuration with uniform entropy per nucleon and chemical com- 
position will satisfy Eq. ([9]) for equilibrium, if and only if the quantity M defined by 



9 

M = Qd-2 J p{r)r^dr is stationary with respect to all variations of p{r) that leave un- 
changed the quantity A^ = VL£,_2 J vAn(r)r^dr and that leave the entropy per nucleon and 
the chemical composition uniform and unchanged. The equilibrium is stable with respect to 
radial oscillations if and only if M or equivalently ii^ is a minimum with respect to all such 
variations. 



These theorems have been proven in general relativity in 4D space-time |28l . |29J . The 
proof of theorems 1 is unconcerned with dimensions and the metric of space-time, therefore 
it is obviously appropriate for our D > 5 star in Gauss-Bonnet gravity. Now we are going to 
prove theorems 2 in present case by using Lagrange multiplier method. M will be stationary 
with respect to all variations that leave A^ fixed if and only if there exists a constant A for 
which M — XN is stationary with respect to all variations. In general, the change in M — AA^ 
for a given variation 6p{r) is 



poo PC 

6M - X5N = VLd-2 / 6p{r)r^-^dr - AfiD-2 / 
Jo Jo 



l + ii_S(r,M) 



- A^._. r nir)SMirydr 

Jo 2{D - 2)anD^2^{r, M)(l + f| - S(r, M)) 2 



where 'E{r^ M) = \/j^ + (£)_2)^n^ r^-s - Here the integrals are carried to infinity for 
notational convenience since they vanish outside a radius R + 6R. These variations are 
supposed not to change the entropy per nucleon, that is = 6{-) + p6{-), therefore 
Sn{r) = p(,"+p(,) '^p(r). Obviously, 6M{r) = ^0-2 Jo Sp{r')r'^^^dr'. Thus, 6M - X6N will 
vanish for all 6p{r) if and only if 

1 _ n 

^ ^ r njr'y^dr' 

"^Jr 2(D-2)(5fi,5-2S(r',M)(l + g-S(r',M))f 

This will be the case for some Lagrange multiplier A if and only if the right-hand side is 



10 



independent of r that 





n' n{p' + p') 1 

^+P~ ^P + P^" ,/l + S-S(r,M) 



+ 



2a 
^ 2a ^ D-2 (D-2)Ci,_2r 



^ D-2 rD-2^Qr. orO-4 rZ^l^r, iWj 



P + P 2aS(r,M)(l + §-S(r,M))^ 



25 
^3 



2iD - 2)aS(r, M)(l + g - S(r, M))l 



(24) 



25 

The condition of uniform entropy per nucleon gives = (^)' + p(^)' and therefore ra'(r) = 
"(r)-f' (1 ■ Substituting the expression of n'(r) in to Eq. (^^, we obtain Eq. (^. Consequently 
SM vanishes for all Sp{r) under a given N and Eq. (jH]). Thus, the stability is determined 
by the sign of second order of 6M. Namely, the star is stable if 5'^M{r)/5p{rY > and 
unstable if 5'^M{r)/5p{rf < 0. 

V. CLOSING REMARKS 

In this paper, we have investigated static spherically symmetric stars in the D > 5 
Gauss-Bonnet gravity without the cosmological constant. It is shown that there is an exact 
solution of such system relating to the exterior solution of the static spherically symmetric 
black hole, and the Gauss-Bonnet interaction decreases or enhances the thermal energy and 
increases the gravitational energy of stars respectively. The Gauss-Bonnet terms will modify 
the upper bound of red shift of spectral lines from the surface of stars of uniform density, 
and the upper bound is dependent on the value of density rather than a constant in general 
relativity case. Besides, two theorems for judging the stability of equilibrium of stars in 
D > 5 Gauss-Bonnet gravity are proved, which are the natural generalization of ones in 
general relativity. 
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